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Semileptonic, B -kW, decays are studied on the MILC dynamical configurations using NRQCD heavy and 
Asqtad light quarks. We work with light valence quark masses ranging between nis and ms/8. Preliminary simple 
linear chiral extrapolations have been carried out for form factors /y and f± at fixed iJ^. The chirally extrapolated 
results for the form factors f+{q^) and fo{q^) are then fit to the Becirevic-Kaidalov (BK) ansatz. Preliminary 
estimates of the CKM matrix element \Vub\ are presented based on recently published branching fractions for 
B" — > 7r~ I'^u exclusive decays by the CLEO collaboration. 



1. Introduction 

First principles calculations of B meson 
semileptonic decay form factors are crucial for 
determining the CKM matrix elements \Vcb\ and 
\Vub\- Recent progress on the lattice towards 
this goal comes from two major developments: 
the ability to go beyond the quenched approx- 
imation with close to realistic dynamical quark 
content [1,2] and the use of improved staggered 
light quarks in heavy- light simulations [3]. We 
report here on unquenched studies of _B ^ tt, IV 
decays on the lattice using one of the coarse 
MILC Nf = 2 + 1 dynamical sets [1], NRQCD 
b quarks and improved staggered (Asqtad) light 
quarks. The light dynamical quark mass is fixed 
at mdyn = 'rns/4: and we vary the light valence 
quark mass between nis and rUs/S. 
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with 

A lattice calculation of the relevant matrix ele- 
ment starts with the three-point correlator 

c(3)(p;,PB,t,rs) = 



2. Form Factors 

Semileptonic form factors parameterise the 
hadronic matrix elements of electroweak currents 
between a B meson and a tt or a p. In particular. 
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where and $b are interpolating operators for 
the TT and B mesons respectively. The three-point 
correlator is fit to the form 
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Figure 1. The form factors f+{q^) and fo{q^), ex- 
trapolated to the physical pion. Only statistical 
errors are shown. The full curves correspond to 
a Becirevic-Kaidalov (BK) parametrization fit to 
the data (see text). 



and the goal is to extract the ground state contri- 
bution Aqq. Working in the B meson rest frame, 
the ground state amplitudes Aoo{V'^) are related 
to the form factors /|| and f± in a simple way. 
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3. The Form Factors /+((7^) and j^{q^) at 
the Physical Pion 

Although our simulations have been carried out 
with light quark masses as low as ms/8 one still 
needs to extrapolate the form factors, determined 
above, to the physical pion. To date, we have 
only carried out simple linear chiral extrapola- 
tions. We first interpolate /|| and /j_ to common 
values of -Ett, the pion energy in the B rest frame. 
These are then extrapolated linearly to the phys- 
ical pion for several fixed values of E^^. 

From the chirally extrapolated /|| and /j^ one 
obtains the form factors f+{q^) and foiq^) at the 
physical pion. This is shown as the circles in 
Fig.l. Our results are currently limited to the 
> 15 GeV^ region. Recently, a very promis- 
ing approach to low form factors has been de- 
veloped, namely "Moving NRQCD" , which will 
allow us to overcome this limitation [5,6]. In 
the mean time, however, we will rely on a model 
ansatz to extend our form factor results into the 
low regime. Specifically, we employ an ansatz 
introduced by Becirevic & Kaidalov (BK) [7], 
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[q = q /Mgt). This ansatz satisfies the kine- 
matic constraint /+(0) = /o(0), HQET scaling 
laws and the requirement of a pole in f+{q^) at 
q^ = Mg,. We find an excellent fit to this BK 
ansatz using the physical Mb' mass and this is 
shown as the full curves in Fig.l (a satisfactory 
BK parametrization was not possible before the 
chiral extrapolation). The fit parameters are 

Cb = 0.42(3), as = 0.41(7), /3s = 1-18(5), 

which translates into 



where Q g are the groundstate amplitudes of the 
TT and B two-point correlators respectively and 
Zvi^ are the matching factors between the lattice 
heavy-light current V^^^ and the current in the 
continuum theory. We use one-loop estimates for 
these matching factors [4]. 



/o(0) = /+(0)= 0.251(15) 

and an effective pole in ,fo{q^) at q^ = (Mjp'*^)^ = 
33.35(1. 36)GeV2. Both /o,+ (0) and MY^" are in 
good agreement with a recent semileptonic B de- 
cay analysis based on Sum Rules [8]. The data 
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points in Fig.l and the /o,+ (0) quoted above in- 
clude only statistical and fitting errors. Further 
systematic errors are discussed in the next sec- 
tion. A comparison of our new dynamical form 
factor results with old quenched data is given, for 
instance, in reference [9]. 

4. Estimating \Vub\ 

The CLEO collaboration has published branch- 
ing fractions for exclusive scinilcptonic B decays, 
including binning into several ranges [10]. We 
combine these experimental inputs with lattice 
results for f+{q^) to extract values for the CKM 
matrix element \Vuh\- The differential decay rate 
for B^^-K-,l+u, 

can be integrated to give ^^-^ and the partial 
width r can be determined from CLEO's branch- 
ing fraction and the Particle Data Group's 
lifetime of 1.542 ± 0.016ps. Our preliminary es- 
timate for I Vub I is then 

r 3.86(32)(58) X 10-3 < < ^2^^^ 
IKfcl = <^ 

[ 3.52(73) (44) x IQ-^ IQGeV'^ < 

where the two values correspond to cither using 
the entire allowed q-^ range or restricting both ex- 
periment and theory to the q"^ > 16 GeV^ region. 
The first error is experimental and the second is 
our current best estimate of lattice statistical and 
systematic errors added in quadrature. In ad- 
dition to 4 6% statistical errors, we estimate 
~ 9% higher order perturbative matching, ~ 5% 
chiral extrapolation, ~ 5% relativistic and dis- 
cretization errors. This adds up to 12.5% lat- 
tice errors for \Vub\ obtained from the q^ > IQ 
GeV^ region. For \Vub\ based on the full q"^ range 
we increase the lattice errors to 15% (an addi- 
tional 8% added in quadrature) taking into ac- 
count the need to rely on the BK parametrization 
to enter the low q"^ region. 

5. Summary 

Unquenched simulations of heavy meson 
semileptonic decays are now feasible and we re- 



port here on the first such calculations using 
NRQCD heavy and Asqtad light quarks (see talk 
by Okamoto for results using Fermilab heavy 
quarks [11]). The use of the improved staggered 
light quark action has allowed for significantly 
smaller statistical and chiral extrapolation errors 
than in the past. Combining lattice results for 
f+iff) with experimental branching fraction data 
has led to preliminary estimates of | Vub \ ■ 
Many improvements are planned: inclusion of all 
dimension 4 (1/M, Us/M and aas) current cor- 
rections, more sophisticated chiral extrapolations 
[12], and simulations at other dynamical quark 
masses and lattice spacings. Use of "Moving 
NRQCD" [5] will also allow us to simulate di- 
rectly at lower q^. 
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